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Abstract 

We will consider a SUSY-S'C/(5) with one right-handed neutrino with a large top like Yukawa coupling. 
Assuming universal soft masses at high scale we compute the low-energy spectrum and subsequently the 
neutralino LSP relic density taking also into consideration SU (5) as well as the see-saw running effects above 
the gauge coupling unification scale. We found that there exists no viable region in parameter space for 
tan/3 $ 35. The f coannihilation process starts becoming efficient for tan/3 ~ 35 — 40. However, this 
process is significantly constrained by the limited range in which the stau is lighter than the neutralino. 
In fact, for a given tan/3 we find that there exists an upper bound on the lightest neutralino mass (A/^o) 
in this region. The A-pole funnel region appears at very large tan/3 ~ 45 — 50, while the focus-point 
region does not make an appearance till large (mp, Mx/2)i namely a few TeV. Large terms at high scale 
can lead to extended regions consistent with WMAP constraints and remove the upper bounds in the stau 
coannihilation regions. 
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I. INTRODUCTION 



The presence of a natural Dark Matter (DM) candidate has been one of the hallmarks of 
low energy supersymmetry models (SUSY). Within the Minimal Supersymmetric Standard Model 
(MSSM), imposition of R-parity conservation would lead the lightest SUSY particle (LSP) to be 
stable and preferably neutral and a color singlet, a perfect candidate to explain the Dark Matter 
(DM) relic density in terms of a WIMP (Weakly Interacting Massive Particle) With the advent 
of precision cosmology, pinnacle of it being the recent results from the WMAP experiment Q], the 
dark matter relic density is now known with high accuracy ^ . Together with availability of state- 
of-the-art numerical tools it is now possible to compute the relic density within a given SUSY 
model accurately up to few percent level. 

However, it has been known that within most SUSY models, the relic density computed is 
either too large or too small; it is required that certain precise or critical relations between the 
various soft SUSY breaking parameters to exist in order the neutralino relic density comes out 
consistent with the measurements. The accuracy with which these relations need to be satisfied 
typically leaves out very tiny special regions in the parameter space where neutralino relic density 
is satisfactory; the tiny regions can be traced back to existence of one of the relations between the 
parameters 4 1 . 

In the simplest models of gravity mediation, i.e, minimal supergravity (mSUGRA), which will 
be the focus of our work, the soft terms have 'universal' boundary conditions at a high energy 
scale, such as the Grand Unification scale, Mgut- The lightest supersymmetric particle (LSP) 
is the lightest neutralino (xi) iii most of the parameter space. There exists three regions where 
viable DM relic density is possible after taking into account the existing low energy direct and 
indirect constraints on the supersymmetric spectrum including the recent constraints from LEP 
measurements These are the (i) The stau coannihilation region, (ii) The A-pole funnel region 
and (iii) The focus-point or the hyperbolic branch region. In the region (i), where the lightest 
stau (fi) mass is very close to the LSP mass 61, significant enhancement in the annihilation cross- 
section happens due to the coannihilation between the lightest neutralino and the fi [7|, |g] (ii) In 
the so-called funnel region, resonant enhancement of the annihilation cross-section takes place in 
the process {xiXi ~^ ff) through the intermediate state of {h, A, H), where h, H, A are represent 
light neutral, heavy neutral and the pseudo-scalar higgs respectively [91]. And finally the focus-point 
region (iii) is the region with large universal soft parameters which however leads to a small weak- 
scale value of ^ (the Higgs bilinear coupling in the superpotential), raising the Higgsino component 
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FIG. 1: Feynman diagrams for LSP ncutralino annihilation into a fcrmion pair (a) and into massive gauge 
bosons and Higgs bosons (b). The relevant parts of the amplitudes are shown explicitly. V and Z are the 
chargino and neutralino mixing matrices. 



in the LSP, and thus leading to significant enhancement in the annihilation cross-section 

A crucial aspect of this arrangement of things where SUSY DM would require some critical 
relations to be satisfied between various parameters is that these relations would not be stable 
under slightest modifications to the SUSY soft masses. These modifications are natural when 
we consider extensions of the minimal supersymmetric standard model (MSSM), where MSSM is 
embedded into either a SUSY Grand Unified Theory (SUSY-GUT) and/or augmented with a see- 
saw mechanism. Additional radiative effects present in these models are sufficient to significantly 
modify the spectrum at the weak scale. More importantly, these effects can destabilize the special 
regions in the parameter space where these relations hold or sometimes these modifications can be 
so large that there might be no parameter space left where these special relations can hold. This 
is precisely the issue we wish to address in our present work. 

We will consider a SUSY-S'C/(5) with one right-handed (RH) neutrino {SU{5)Rj\f) with a large 
top like Yukawa coupling^. Such a model is motivated from the requirement of non-zero neutrino 
masses and can have UV completion in 50(10) based models. Assuming universal soft masses at a 
scale Mx — 10^^ GeV, we compute the low-energy spectrum taking also into consideration SU{5) as 
well as the see-saw running effects above the gauge coupling unification scale ( Mqut — 10^^ GeV). 

^ We need at least one more RH neutrino and a corresponding Yukawa coupling to fit the neutrino data. This 
Yukawa coupling could be much smaller than the top Yukawa and thus its effects negligible to the discussion 
pertaining here. 
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Adding the effects of the further running, down to the weak scale, we compute the rehc density 
with the weak scale spectrum. 

Firstly, we found that there exists no viable region in parameter space for tan/? ~ 35 for 
neutralino dark matter. Parameter space points consistent with WMAP results start appearing 
with tan/3 ~ 35 and increase with increasing tan /3. However, this process is highly constrained 
by the limited range in which the stau is lighter than the neutralino. In fact, we find that the 
coannihilation region is significantly modified with respect to the MSSM and it has the shape of a 
trunk, in the {mQ, Mij^) plane^, cutting into the M112 axis. Thus, for a given tan/3, there exists 
an upper bound on the lightest neutralino mass (M^o). The allowed parameter space widens up 
at very large tan /? ~ 50 where the A-pole funnel region appears. The focus-point region is very 
difficult to realise in this model because electroweak symmetry breaking is much easier to achieve 
compared to the standard MSSM. In fact, even going up to 5 TeV in (mo,M]^/2) plane, focus-point 
region does not make an appearance. 

While the above results are particularly true as long as the universal trilinear coupling Aq = 
at the high scale, the results can significantly change with large —3 mo ^ Aq ^ 3 mo boundary 
conditions which due to the possibility that some of the masses becoming tachyonic, significantly 
affects the viable parameter space. For instance, for Aq = +3 mo, the stau coannihilation region 
branches out into two regions, a new region other than the trunk region discussed above. The 
A-pole funnel also spreads giving rise to a larger region of the parameter space which is viable 
under WMAP constraints. 

The rest of the paper is organized as follows: in the next section, we will elaborate our framework 
along with revising some aspects of neutralino dark matter. In section 3, we summarize the 
constraints we use in our numerical analysis as well as sketch our computing procedure. In section 
4, we discuss our results and close with a discussion in our final section. 



II. SU{5}rn and NEUTRALINO DM 

As mentioned in the introduction, in the following we will consider a simple SUSY-GUT frame- 
work based on SU{5). Further we also add a RH neutrino with a large top-like Yukawa coupling. 



which is natural when this model is incorporated in an 50(10) model 



The superpotential of 



mo and Mij2 are the universal mSUGRA boundary conditions for the scalars and gauginos masses, respectively. 
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this framework is given by: 
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(1) 



i,j being generation indices and 10 and 5 (5), representing the tensorial and fundamental (anti- 
fundamental) representations of SU{5) respectively. The subscript u, d denote the up and down 
type Higgs fields^. Two comments are in order before we proceed further. 

Firstly, to accommodate neutrino masses, we would need at least one more RH neutrino in 
addition to the one already present in Eq. ([1]). However, in most models, the corresponding Yukawa 
couplings are not as large as the top Yukawa coupling^. Thus, the impact of the second and third 
neutrino Yukawa couplings even if they are present is expected to be minimal on the discussion 
below. In fact, we have explicitly checked this to be the case in models where the neutrino Yukawa 
coincides with the up-type Yukawa matrix. Secondly, keeping 5'O(10) completion in mind, we have 
chosen the high scale Mx to be around ~ 10^^ GeV. Another natural choice would be the reduced 
Planck scale, where Mx ~ 10"'^'^ GeV. Nevertheless, with this choice of the scale, we believe our 
limits on tan/3 (as well as the neutralino mass) will only become stronger. In Sec. llVl we will briefly 
comment about the dependence of our results on the choice of Mx- For most of the discussion, 
however, we will just choose the lower value for Mx- 

The soft SUSY-breaking terms are considered to be universal at the scale, Mx- These contain 
the standard mass-squared terms for the scalar particles (the squarks and the sleptons), the bilinear 
Higgs scalar couplings, the trilinear scalar couplings as well as the mass terms for the gauginos. 
Given that R-parity is conserved, the LSP remains stable. In most models of SUSY-breaking like 
mSUGRA, the LSP is the lightest neutralino. The neutralinos are the physical superpositions of 
the fermionic partners of the neutral electroweak gauge bosons, called bino {B^) and wino iW^), 
and of the fermionic partners of the neutral Higgs bosons, called Higgsinos (-ff^, H^). We can 
express the lightest neutralino as 



where Zij are the elements of the neutralino mixing matrix. It is commonly defined that Xi is 
mostly gaugino-like \i P = \Zii\'^ + \'^\2^ > 0.9, Higgsino-like if P < 0.1, and mixed otherwise. 
This factor P has a major impact in determining the relic abundance of the LSP as it determines 
the rate of neutralino annihilation cross-section. In Fig. [1] we show the relevant Feynman diagrams 



X? = ZiiB + Z12W3 + ZnHll + ZuH\ 



u ' 







(2) 



^ We have written the superpotential in the basis where up-type quarks are diagonal. 

* In fact, if one considers models based on 50(10), it can be as small as the charm Yukawa, j/c- 
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FIG. 2: Schematic picture of the energy scales involved in the model. 

contributing to neutralino annihilation. This can be computed for a given SUSY spectrum, which 
depends on the SUSY breaking model. 

In the MSSM, with universal soft terms at Mgut = 2 x 10^^ GeV (CMSSM), for most of the 
parameter space, the lightest neutralino is mainly bino and, as a consequence, the annihilation 
cross-section is small producing too large relic abundance. Nevertheless, as mentioned in the 
introduction, three corridors exist where the cross-section is enhanced and the WMAP bounds can 
be satisfied. First, there is the coannihilation branch, i.e. the region where the stau mass is almost 
degenerate with that of LSP. Another corridor corresponds to the A-pole region. This occurs in 
the parameter space where 4(m^o)^ — ~ ''^Hi ~ "^^2 ~ -^1' '^'^^^ being the soft Higgs 
mass terms, and the dominant neutralino annihilation process is through the 'resonant' s-channel 
pseudo-scalar Higgs exchange. Finally, there is a 'Higgsino' corridor close to the no electroweak 
symmetry breaking (No EWSB) region. In this region the fi parameter is much smaller than the 
bino mass and the lightest neutralino has a larger Higgsino component (thus annihilating efficiently 
via Z boson exchange). 

As we have stressed in the introduction, these corridors or special regions in the parameter space 
are not stable under modifications to the soft spectrum. In this connection, it was already noted, 
for e.g. in Refs. [1^, that the annihilation cross-section can change significantly if the soft terms 



13l | assumed no-scale 



determined at the unification scale, MquT; are non-universal. Similarly, Ref. 
boundary conditions at Mx = Afpiancki which would result in non-universal boundary conditions 
at Mgut^- In the scheme we are considering, which has been elaborated in the beginning of this 
section, something similar happens: the soft terms are indeed non-universal at the GUT scale and, 
as we will discuss in detail below, this would have implications for the regions in the parameter 
space where viable relic neutralino density was possible. 



Our work is closest in spirit to this work and could be considered a generalization of it. 
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III. SPECTRUM AND CONSTRAINTS 



Before proceeding further in the present section, we will provide details of our computation of 
the supersymmetric spectrum in the framework we have considered and the various constraints 
we have used. A schematic diagram describing the integration procedure is presented in Fig. [2l 
As inputs at the weak scale, we consider the Yukawa couplings of the up~type quarks, down-type 
quarks, charged leptons and tan /3. We consider one generation of neutrino with a mass around 0.05 
eV. We use 1-loop RGEs to run all the Yukawa couplings up to the high scale. A more detailed 
description of the procedure can be found in Ref. [l^. After several iterations to check consistencies 
between neutrino-scale spectrum and Mx-, we are now ready to compute the spectrum from the 
high scale to the weak scale. The soft parameters, gauge and Yukawa couplings follow the SU(5)j:j^ 



RGEs [15|] between Mx and Mqut- We then compute the evolution of these SUSY parameters 
using a modified version of the Fortran package SUSPECT 1^, where we have included one RH 
sneutrino between the GUT scale and the RH neutrino mass scale, Mr. The runnings then follows 
the classical MSSM RGEs after Mr where the RH neutrino sector decouples. The low energy mass 
spectrum is then calculated using SUSPECT. The evaluation of the b ^ branching ratio, the 
anomalous moment of the muon and the relic neutralino density is carried out using the routines 
provided by the program micrOMEGAsl .3.1 



18l ]. The following are the list of constraints we apply 



on the soft spectrum: 



(i) Electroweak symmetry breaking : 

Minimizing the Higgs potential in the MSSM leads to the standard (tree-level) relation 

„ — m?r tan^/3 + m?r 1 „ 

This minimization condition is imposed at the scale MsusY = yj'^ii "^fe ' "^^^^re m^, are the stop 
masses. Eq. ([3]) can be approximated in most cases by 

-ml^ - (4) 

When the right-hand side is negative, electroweak breaking cannot occur. The Higgs mass pa- 
rameter rrijj^ is positive at the GUT scale, but decreases with decreasing scale down to MsusY, 
through the contributions it receives from RG running at the scale jl; in presence of RH neutrinos, 
the RG equation for mj^^ reads: 

^^''^ %in(/2/Mx) ^""^2 + + "'Oa + + '^vK^k + ^% + ^l, + (5) 



where gauge contributions have been neglected. Here and m% corresponds to the ii and 
ifi soft masses, while m? , m?- are the same for the left-handed and right-handed sneutrinos 

' L3' Af ° 

respectively; yt and yi, are the top and the neutrino Yukawa couplings, while At and A^, the 
corresponding soft SUSY-breaking trilinear couplings. Typically, in the MSSM the value of m'j^^ 
at the scale MgusY depends mainly on the gluino mass M-^ at the GUT scale (via its effect of 

increasing m% and rn^ ). However, in our case, the further increasing of m% and rn^ due 

Us Qs U-i Qz 

to the GUT running between Mx and AfcuTi together with the positive additional contribution 
proportional to y^, drives rn?i{^ toward even more negative values. As a consequence, the parameter 
space region where /i^ < (or No EWSB region) is much more reduced, and even absent in our 
case as we will see in the analysis. 



(ii) The mass spectrum constraints: 

We have implemented in our analysis the lower bounds on the masses of SUSY particles and 
of the lightest Higgs boson. In the squark and slepton sector we checked for the occurrence of 
tachyons. We applied in our analysis the LEP2 lower bound limit on the mass of the lightest 
chargino m-+ > 103.5 GeV. In the non-tachyonic region, typically, the most constraining is the 
lightest Higgs boson mass constraint. In the decoupling limit {Ma ^ M^, applicable in all our 
parameter space), m/^ > 114 GeV at 3cj. This bound is very sensitive to the value of the top mass. 
We have taken mt = 174 GeV throughout our analysis. 



(Hi) The b ^ s^ branching ratio: 

One observable where SUSY particle contributions might be large is the radiative flavor changing 
decay b — > 57. In the Standard Model this decay is mediated by loops containing the charge 
2/3 quarks and VF— bosons. In SUSY theories additional contributions come from loops involving 
charginos and stops, top quarks and charged higgs bosons or, in presence of RG induced mixing 
in the down squarks sector, down squarks and gluinos. The measurements of the inclusive decay 
B at CLEO Q and BELLE Q, 

leads to restrictive bounds on the branching ratio 
b — > 57. We impose in our analysis 2.33 x lO"'^ < BR{b 57) < 4.15 x 10~^ at the 3fT level. 
We mostly choose fi > enabling cancellations between chargino contributions and charged Higgs 
contributions [21 1. 



Before proceeding further, we elaborate a bit more about the flavour processes in our model. 
Since we are in a Grand Unified theory, one would expect large off-diagonal entries to be generated 
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in the soft matrices at the GUT scale due to RG running effects from Mx to Mqut- This could 
have large effects in flavour processes at weak scale. For example, a large off-diagonal entry in the 
23 sector of the down squarks would contribute dominantly through the gluino interactions to the 
6 — > s + 7 rate leading to significant constraints on that entry, unless of course, there are accidental 
destructive interferences among the various contributions already mentioned above. However, in 
our case, even if the corresponding mass matrices are non-universal at the GUT scale, the entries 
are generically small as they are proportional to the CKM entries. Some ambiguity can arise in 
the way we handled the neutrino sector. If all the three right handed neutrinos are present and the 
neutrino Dirac Yukawa matrix carries large "left" -mixing then they can generate large off-diagonal 
entries at the Mji scale. In such a case a combination of both leptonic and hadronic flavour violating 



constraints have to be used 



22l |. On the other hand, it is known that even small departures from 



universality at high e nerg y may have a strong effect on 6 — > s + 7, especially tending to make 



the constraint weaker 



23( 1 . In the present case, the constraint will be a bit weaker because RH 



neutrinos in SU{5) only affect the down squarks 'RR' sector. Moreover, we can avoid these effects 
at least in our framework by choosing small left mixing, though flavour violation in the leptonic 
sector might still be significant even in that case [l^. [2^. We postpone such an analysis for future. 



(iv) The anomalous moment of the muon: 

We have also taken into account the SUSY contributions to the anomalous magnetic moment of 
the muon, = {g^ — 2)/2. We used in our analysis the recent experimental results for the muon 



anomalous magnetic moment 
Model contributions 



25( 1 , as well as the most recent theoretical evaluations of the Standard 



26| . It is found that when e'^e data are used the experimental excess in 



{g^ - 2) would constrain a possible SUSY contribution to be 7.1 x 10~^° ^ af^^'^ ^ 47.1 x 10"^° at 
2(7 level. However when tau data is used a smaller discrepancy with the experimental measurement 
is found. In order not to exclude the latter possibility, when analyzing the parameter space with 
^ > we will simply plot contours with the relevant value a^^^^ = 7.1 x 10"^'^. 



(v) Relic Density: 

Our basic assumption is that the LSP is stable on cosmological time scales. Furthermore we will 
assume that the LSP abundance is thermal. Within such framework the regions of the parameter 
space that lead to overproduction of dark matter are excluded. On the other hand, the regions 
that yield LSP abundance below the WMAP limit are not considered as excluded (though as less 
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FIG. 3: The {mo, M1/2) plane with aU the low-energy constraints on the parameter spaee. While the color 
code is explained above, the green dashed line indicates the Higgs mass bound from the LEP, while the dark 
dashed line represents the 6 — *■ s + 7 limit. 



favored), but simply require non-thermal production or a dark matter candidate beyond the soft 
spectrum. The WMAP collaboration gives the 3a narrow limit 

0.087 < ^DMh^ < 0.138 (6) 

on the dark matter relic abundance. 



IV. SU{5)rn and CiDM 

Before going to present the numerical results of our analysis, let us briefly mention the major 
changes to the spectrum one could expect in our framework and the possible implications on the 
neutralino relic density. As we have mentioned in the previous section, we have assumed that the 
soft masses are universal at a scale Mx which we have chosen to be around 10"'^'^ GeV. As far as 
the neutralino relic density is concerned, two major effects on the soft spectrum due to the SU{5) 
running between the scale Mx and Mqut and subsequently to Mr can be given as [l^ : 

• The right-handed slepton fji now sits in the 10 and thus it will receive contributions from 
the full gaugino multiplet of the SU{5). At the leading log level, these contributions are 
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given at the GUT scale as : 

mljMGVT) ^ ^«5M,%ln(^^) ^ 0.25 Af,%, (7) 

where we have taken the Hmit tuq — > and as = |^ represents the unified gauge coupfing 
at Mgut- This large positive contribution from to the tr make the stau heavier than the 
LSP (which is the lightest neutralino) for most of the parameter space. Let's note here that 
choosing a larger value for Mx (such as Mpianck) would further increase rrif^, strengthening 
the effect described above. 

• In our framework the neutrino Yukawa coupling taking part in the see-saw mechanism as 
shown in Eq. ([T]) is taken to be as large as the top Yukawa coupling. This introduces an 
additional top- Yukawa like coupling to the up- type Higgs from the scales Mx down to M/j. 
Compared to the case of MSSM with RH neutrinos, as shown in Eq. ([5]), here the length 
of running is more and thus more effective radiative electroweak symmetry breaking takes 
place^. Moreover, the increasing of the up-squarks soft masses, due to the unified gauge 
sector, contributes to pushing m'jj^ down to more negative values than in the CMSSM. As a 
consequence, for most of the parameter space, we find that electroweak symmetry breaking 
is viable unlike in the CMSSM where substantial amount of parameter space is ruled out. 

The above two reasons are sufficient to offset the conditions which give viable dark matter in 
CMSSM at low tan /3. As stated in the introduction, we found no viable region giving the required 
relic density up to tan /? ~ 35. In Fig. [3l we show the non-existence of any region of the parameter 
space which gives rise to the correct relic density in our framework of SU (5) /jjvi for tan /? = 10 and 
^0 = 0. All the three branches of CMSSM are not possible here^ as none of their corresponding 
relations could be satisfied within this regime. 

To make this statement more concrete, in Fig.^J we plot all the points which satisfy the available 
direct /indirect constraints and give viable relic density as a function of tan/3 and LSP mass. Prom 
the figure we see that viable DM is only possible for values of: 

tan/3 ^ 34 ; m^o ~ 150 GeV (8) 

^ The value of Mr is uniquely fixed {Mr ~ 6 x fO^^GeV), through the see-saw mechanism, by our choice of the 
neutrino Yukawa coupling and light neutrino mass scale. Some variation in AIr can come in a complete three 
generation model choosing inverse/normal hierarchy and particular mixing patterns. Choosing a larger light 
neutrino mass {0{eV)) would simply lower Mr, thus increasing the RG effect on ji. 

^ Even if we have confined ourselves, in the figure, to the parameter space below f TeV for both mo and A/1/2, we 
have performed a scanning up to 5 TeV, without finding viable regions. 
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FIG. 4: Points allowed by experimental and theoretical constrains after a scan on (0 < mo, M1/2 < 1 TeV) 
and (20 < tan/3 < 55) 

These are quite strong lower bounds on the neutralino mass and tan/3 and will be useful in dis- 
tinguishing this model compared to the standard CMSSM parameter space. In the Fig. [H we see 
that the small strip of points on the lower left corner gives a more constrained solution in m^o 
for tan/? ~ 45. For larger tan/?, we see that almost all possible values are allowed in the lightest 
neutralino mass. However, at this level, we are unable to distinguish between the actual annihila- 
tion processes involved in generating the correct relic density. A more detailed look shows us the 
existence of two regions which we elaborate here. 

A. Coannihilation Region 

As we have discussed previously, the large pre-GUT scale contributions to the in Eq. ([7]) 
make the lightest stau heavier than the LSP in most of the parameter space. However, when tan /? 
becomes large, it is possible to make the lightest stau fi closer to the LSP mass by making the 
left-right (LR) mixing term mT-{Ar — fitanP) large. This is what indeed happens in our case. 
Stau coannihilation cross-sections start becoming effective for tan/3 ^ 34, leading to relic density 
within the WMAP bound. In Fig. [5] we have plotted, the WMAP compatible regions for the case 
of CMSSM as well as SU{5)rn, for tan /? = 40 and Aq = 0. 

Two features are evident from such a comparison between the CMSSM and the SU{5)rn- (i) 
the WMAP compatible region is much smaller in the case of SU{5)rn and (ii) the shape of the 
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FIG. 5: Wc show the coannihilation region for tan/3— 40 for the case of CMSSM and SU{5)rn- The green 
dashed and the dark dashed hnes represent same as above. 

allowed region is quite different in the case of SU{5)rn (we will call it trunk region). In fact, the 
allowed region cuts-off for a value of M1/2 — 520 GeV in the above plot. This would correspond 
to a LSP mass of around 240 GeV; an LSP of higher mass would not be able to give rise to the 
correct stau-neutralino coannihilation rate. This we would say corresponds to an upper bound on 
the LSP mass which could give viable stau coannihilation rate such that the relic density is within 
the WMAP bound. The presence of such upper bound in the LSP mass is not peculiar of the model 
we are considering here, since similar bounds are well-known also within CMSSM [3, Is^]- However, 
the nature of these upper bounds is quite different in the two models. In the CMSSM, the reason 
is that the coannihilation cross-section decreases at larger SUSY masses [8|] and the coannihilation 
rate becomes too low even if the condition m^^ ~ m^o is satisfied. On the other hand, in SU{5)rn, 
the region for which m^^ ~ m^o simply disappears for some value of M1/2, as a consequence of the 
peculiar shape of the f-LSP region. This excludes the possibility of efficient coannihilation even 
for values of tuq and outside the "1 TeV box" of Fig. [5] (as we numerically checked), so that 
the value m^o ~ 240 GeV is a real upper bound for the case considered (tan P = 40, ^0 = 0). 

The reason for these peculiar features of SU {5)^^ can again be traced back to the large gaugino 
contribution to the stau mass above the GUT scale. At the weak scale, roughly the stau mass matrix 
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is now given by: 

Mi = I | , (9) 

where, including the pre-GUT effects, 

m?^^ ^ (l-;9)m2 + 0.3M2/2 (10) 

with p being a positive coefficient dependent on tan/3 and smaller than one in the present case^, 
and 

m|^^ = rrir {A^- — p.ta.n fi) 

~ —irij- fitanP (11) 

where we explicitly indicate the dominance of the tan /3 enhanced term. In first approximation the 
lightest eigenvalue of Eq. ([9]) is given by: 

~ THf^^ — rrir p tan P 

~ {1- p)ml + 0.3 - mr 1^ tan /3 (12) 

The coannihilation condition requires that niri should be almost degenerate with the mass of the 
LSP, which is mostly the Bino mass. Mi ~ 0A7 A-Ii/2- Taking mo ~ -^^1/2 ~ ~ -^^SUSY and 
barring the tan/? dependence of the coefficient p, we find as a result tan/3 ~ (1 — p){Ms\jsY /t^t), 
from which, for reasonable values of MsusY ~ 200 GeV and considering that in the present case 
0.1 :$ p ^ 0.7, we get the lower limit tan/3 ~ 35^. Solving numerically the system with the 
parameters varying in a wide region of the parameter space (0 < niQ, M1/2 < 1 TeV, 5 < tan /3 < 50, 
Aq = 0), we find that that there is no solution for tan/3 ~ 27, furthermore, for a given tan/3, only 
a restricted region of the parameter space in M1/2 is allowed. In Fig. SI this region corresponds to 
the lower strip of points 34 $ tan /3 $ 45. As we can see, the limit on tan /3 is greater after applying 
the experimental constrains mainly due to restrictions coming from the bound on nih- Overall, the 
spectrum prefers a neutralino mass of at least of 0(150 GeV) and further with an upper bound of 
about 0(200 GeV) if stau coannihilation is the dominant process in efficiently reducing the relic 
density to the observed levels^''. 



* On the contrary, when \Ao\ 7^ the parameter p could be driven larger than one by running effects due to the 
A-term. This could cause the arising of tachyonic stau masses, as we will discuss in Sec. II V Dl 

^ A longer running, as in the case Mx = Afpianck, would increase rni^^ in Eq. (|12|l . As a consequence an even 
larger value of tan j3 would be necessary to satisfy to coannihilation condition (though the larger m^^^ would be 
partially compensated by the simultaneous enhancing of Mi). 

Such a narrow range can perhaps be useful to distinguish this model compared to standard CMSSM at colliders. 
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FIG. 6: Same as above for tan/? = 50. 

B. A-pole funnel 

Let us now move to the funnel region. As in the CMSSM, this region makes its appearance for 
large tan /3 ~ 45 — 50. In Fig. [6l we have plotted the funnel region for tan [3 = 50, Aq = 0, again 
in comparison with the CMSSM case. We can also see large regions of the parameter space where 
the lightest stau is the LSP and thus there are regions of coannihilation also which are fused with 
the funnel region. The trunk region is no longer present and neither do the upper bounds. 

C. Focus Point region 

As discussed in Sec. IIIH the REWSB is very easy to achieve due to the presence of additional 
GUT effects and top like Yukawa contribution from the neutrino Yukawa coupling to the up type 
Higgs between Mx and Mr. The result is there is no focus-point region present. We could not 
find the focus-point until (mo,Mi/2) — 5TeV. We have not explored further, since these regions 
would be far from the reach of LHC. 

D. I^ol - 3 mo 

All the solutions discussed above are for the special case where = at the high scale. Given 
the constrained solutions we have here specifically for the stau coannihilation process, which can 
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depend significantly on the non-zero Aq parameter, it is instructive to consider the extreme cases 
of Aq = —3 mo and Aq = +3 itlq. For = ~3 tjiq, we found that the occurrence of large regions 
of the (m-o, M1/2) plane excluded by tachyonic makes such case really disfavoured, leaving 
only a small portion of the parameter space viable. On the other hand, for Aq = +3 tjiq such 
excluded region results much smaller and localized in the portion of the plane where mo is large 
and small. Such behaviour can be explained by looking again at Eq. (jl2p . If the parameter 

p there, which depends on and A^, becomes greater than one because the additional A-term 
running effect with respect to the case ^0 = 0, the coefficient in front of uiq becomes negative. 
Therefore, values of uiq larger than can easily make the lightest slepton tachyonic. The effect 
is enhanced by growing tan /3 both because the direct tan P dependence in Eq. (jl2p and because the 
tan /3-enhanced and contribution to the value of p. The different behaviours for Aq = ±3 mo 
are due to the renormalization group running of A^ itself, which is not independent of the sign of 
Aq. In Fig. [71 we plot, for the Aq = +3 mo boundary condition, the regions which have the relic 
density consistent with WMAP measurements. We see from the left figure that Aq = +3 mo is 
very seriously constrained in the (tan/3, M^o) plane compared to the ^o = case. The appearance 
of an upper bound (~ 45) on tan/? is due to the increasing of the mentioned tachyonic-f region 
for large values of tan/?, as it appears evident from Eq. (jl2p . Moreover, a closer look reveals a 
different branch of the coannihilation region in addition to the trunk region. This is shown in the 
right panel of Fig. [7l We see in addition to the standard trunk branch, a thin strip on the top left 
corner, which appears along the f-LSP branch not present for ^o = 0. Such branch runs along the 
large tachyonic-f region which now appears due to the non-zero A-term as explained above. 



V. REMARKS 



If supersymmetric standard models do really provide a dark matter candidate whose relic density 
constitutes most of the energy budget of the universe, one would typically expect non-trivial 
relations to hold true between various supersymmetric parameters. The point we tried to address 
was how stable these relations are under modifications to SUSY standard models. SUSY-GUTs 
with see-saw are natural extensions of MSSM in the light of experimental evidence for non-zero 
neutrino masses. In the simplest model we have considered here, SU{5)jiN, we found that relic 
density requirements put severe constraints on the range of tan /?. Only large tan /? ~ 35 is now 



allowed. Such large tan/3 will also lead to large effects in flavour physics 
present in this class of models. 



2j| which will also be 
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FIG. 7: j4o = 3 mo case. Left panel: points allowed by experimental and theoretical constrains after a scan 
on (0 < mo, Ml/2 < 1 TeV) and (20 < tan/3 < 55). Right panel: the (mo,Mi/2) plane plot for tan/3 = 40; 
we can see the two branches of the stau coannihilation region. 

Of the three possible corridors of viable parameter space present in CMSSM only two appear 
here. The stau coannihilation region is both bounded from above and below as long as Aq = 0. 
The A-pole funnel appears for large tan j3, whereas focus-point region does not make an appearance 
in this class of models. 

In the recent times, there has been a large interest in tr ying to reconstruct the dark matter 
relic density from SUSY spectrum measurements at colliders 



291 ] ■ While LHC alone might not 



be able to reconstruct the spectrum, one expects that astrophysical relic density measurements 
and LHC data would give a rough handle on the lightest neutralino mass. The couplings and 
cross-sections leading to precision measurements of the SUSY soft spectrum are possible at the 
ILC (International Linear Collider) operational in two modes with -^i = 500 GeV and 1 TeV. In 
particular, the use of r polarization [30(] should be able to distinguish this model from CMSSM 
as long as SUSY spectrum lies in the stau coannihilation region [3l|. Here we use the fact that 
in our model, the relative sizes of m? , m? and m? change with respect to the CMSSM 
and consequently the large left-right mixing of stau required to get the degeneracy with the LSP. 
Further, it might be possible that the trunk region of this model, which has upper and lower 
bounds in the neutralino mass, can be tested from the data at LHC. It is not clear if the other 
regions which can have degenerate solutions with CMSSM can be discriminated either at LHC or 
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a combination of LHC and future linear colliders. Such a question needs further studies. 

Let us conclude by stressing that the peculiar phenomenology of the coannihilation we presented 
here turns out to be a quite general feature of GUT models. In fact, the crucial point is the 
increasing of the lightest stau mass due to effects of the unification gauge group such as in Eq. ([7|) . 
Similar effects are present every time m?^ feels the presence of a unified gauge sector, independently 
of particulars of the GUT structure and of the presence of right-handed neutrinos as in our case. 
Therefore, we claim that even in the absence of see-saw or large neutrino Yukawa couplings, the 
phenomenology of stau coannihilation regions would be similar in most GUT models. The presence 
of large see-saw couplings would effect the radiative electroweak symmetry breaking part and thus 
would have implications for the focus point regions. In the end our work shows that dark matter 
seems to predict large tan/3 in SUSY-GUTs. 
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